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Abstract 

We study the dynamics of an electron wave-packet in a two-dimensional square lat- 
tice with an aperiodic site potential in the presence of an external uniform electric 
field. The aperiodicity is described by em = F cos (vram^"' ) cos (vram^^ ) at lattice 
sites {mx,my), with ira being a rational number, and Vx and i^y tunable parame- 
ters, controlling the aperiodicity. Using an exact diagonalization procedure and a 
finite-size scaling analysis, we show that in the weakly aperiodic regime {vx, < 1); 
a phase of extended states emerges in the center of the band at zero field giving 
support to a macroscopic conductivity in the thermodynamic limit. Turning on 
the field gives rise to Bloch oscillations of the electron wave-packet. The spectral 
density of these oscillations may display a double peak structure signaling the spa- 
tial anisotropy of the potential landscape. The frequency of the oscillations can be 
understood using a semi-classical approach. 
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1 Introduction 

Materials with restricted geometry, such as semiconductor quantum-well struc- 
tures, [1] quantum dots and wires, [2,3] organic thin films [4] as well as 
quasiperiodic structures, [5] are nowadays subjects of growing interest from 
both fundamental and practical points of view. An attributive peculiarity of 
almost all of them is the presence of disorder, which can be both of an intrinsic 
nature (imperfections of the structure itself) and originated from a random 
environment. 

Whenever disorder is involved, Anderson's ideas about localization of quasi- 
particles states come into play. [6] In three dimensions (3D) , the states at the 
center of the quasiparticle band remain extended for a relatively weak disor- 
der (of magnitude smaller than the bandwidth), while the other states (in the 
neighborhood of the band edges) turn out to be exponentially localized. This 
implies the existence of two mobility edges which separates the phases of ex- 
tended and localized states. [7] On the contrary, uncorrelated disorder of any 
magnitude causes localization of all one-particle eigenstates in one dimension 
(ID) [8] and two dimensions (2D). [9] 

Since late eighties, however, it has been realized that extended states may sur- 
vive on ID systems if the disorder distribution is correlated. [10,11,12,13,14,17,15,16,18] 
Thus, a short-range correlated disorder was found to stabilize the extended 
states at special resonance energies. In the thermodynamic limit, such ex- 
tended states form a set of null measure in the density of states, [10,11,12,13,14] 
implying the absence of mobility edges in these systems. In contrast, systems 
with long-range correlations of disorder support a set of delocalized states 
within a finite bandwidth, [15,16] giving rise to mobility edges. Theoretical 
predictions of localization suppression on ID geometries, due to correlations 
of the disorder distribution, were confirmed experimentally in semiconductor 
superlattices with intentional correlated disorder, [17] as well as in single-mode 
waveguides with correlated scatterers. [18] 

Among ID models with extended states, aperiodic Anderson models [19] with 
an incommensurate site potential represent a class of particular interest. These 
models have been extensively investigated in the literature, [20,21,22,19,23] 
and the localized or extended nature of the eigenstates has been related to 
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general characteristics of aperiodic site-energy distributions. Interesting real- 
izations of this kind of distribution are achieved when weakly aperiodic poten- 
tials are considered: these systems support the occurrence of mobility edges, 
in a close analogy with the standard 3D Anderson model. [19] 

A key question in any general theory of aperiodic systems is the relationship 
between a given aperiodic site-energy distribution and physical properties of 
the system. Nowadays, the solution of this problem is still lacking. [24] Never- 
theless, some important results have been obtained concerning the dynamics 
of elementary excitations, like electrons, phonons, excitons, polaritons, spin 
waves, plasmons or magnons, propagating through different classes of aperi- 
odic lattices. [25,5,26,27] In these studies, relevant physical properties were 
analyzed in terms of appropriate model Hamiltonians. 

Recently, the dynamics of a single electron in homogeneous systems subjected 
to a uniform static electric field has received much attention. Under this con- 
straint, the electron undergoes the so-called Bloch oscillations, [28,29,30,31,32] 
whose amplitude is proportional to the energy bandwidth. Electronic Bloch 
oscillations were observed for the first time in semiconductor superlattices [33] 
(see Ref. [34] for an overview). A similar phenomenon of sustained oscillations 
of the electromagnetic field, named photon Bloch oscillations, have been also 
found in 2D waveguide arrays and optical superlattices based on porous sil- 
icon. [35] In disordered systems, Bloch oscillations were predicted to occur 
in both ID [36] and 2D [37] geometries whenever the disorder distribution is 
long-range correlated. The amplitude of the oscillations was found to carry in- 
formation about the energy difference between the two mobility edges, existing 
in these systems. The occurrence of Bloch oscillations was also theoretically 
argued for ID chains with aperiodic slowly varying potentials. [38] A spectro- 
scopic manifestation of Bloch oscillations is found in the appearance of the 
Wannier-Stark ladder in the absorption spectrum. [39] 

In this paper, apart from the generalization of the previous work, [19] done for 
a one-dimensional system, to a 2D system, we consider an electron subjected 
to an external uniform electric field. The latter brings new features to the 
electron wave packet dynamics moving in incommensurate potential. In par- 
ticular, despite a quasi-stochasticy of the system, sustained Bloch oscillations 
can be observed, previously believed to exist only for a periodic structures. 
In order to produce an aperiodic (2D) site potential, the formalism used in 
Ref. [19] was considered. It consists in using a sinusoidal function whose phases 
and (f)y varies as a power-law. The power-law exponent controls the degree 
of aperiodicity in the site potential. We start by employing an exact diagonal- 
ization formalism to compute the dc conductance and participation number 
in the absence of the electric field. We numerically demonstrate that their size 
behavior is consistent with the existence of extended states near the band cen- 
ter for weakly aperiodic potentials. In this regime, the wave-packet dynamics 
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becomes ballistic. Furthermore, we focus on the wave-packet dynamics in the 
presence of an external uniform electric field. We show that in the limit of 
a weakly aperiodic site potential, the electric field promotes sustained Bloch 
oscillations of the electronic wave-packet. The frequency of these oscillations 
can be understood within the framework of a semi-classical approach. 

The outline of the paper is as follows. In the next section we present our model 
and quantities of interest. In Sec. 3, the results of the numerical simulations 
of the dc conductance, participation number and electronic wave-packet dy- 
namics are discussed. We summarize in Sec. 4. 



2 Model and formalism 

We consider a tight-binding single electron Hamiltonian on a regular 2D open 
lattice of spacing a with an aperiodic site potential and a uniform static electric 
field [32,40] 

H^J2{^m+U ■ m) \m){m\ 




where \m) is a Wannier state localized at site m = rrixex + rriyey, is its 
energy, and U = eFa is the energetic bias given by the electric field F, —e 
being the electron charge. Here and Sy are the corresponding Cartesian unit 
vectors. We will assume that the electric field F is apphed along the diagonal 
of the square lattice, i.e., U — U{ex + ey)/\^, were U is the bias. Transfer 
integrals are restricted to the nearest-neighbor interaction and are given by J. 
Hereafter, the energy scale is fixed by setting J = 1. 

The 2D aperiodic potential Cm is taken in the form: 

Em — V cos ijiarny^) cos {Tiarriy^) , (2) 
where V , a, v^i and Vy are variable parameters. 

Equation (2) is a 2D generahzation of the aperiodic ID potential introduced 
by Das Sarma et al. [19] In ID and at = 1, this potential represents just 
Harper's model. In this case, a rational a describes a crystalline solid, whereas 
an irrational a yields an incommensurate potential. It was shown in Ref. [21] 
that for z/ > 1 all one-electron states are localized. In this regime, the frequency 
of the potential oscillations diverge in the thermodynamic limit, resulting in a 
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fairly uncorrelated sequence of the on-site energies. Because of that, this case is 
usually named as the pseudo-random regim,e. Oppositely, it was demonstrated 
that the range Q < v < 1 supports a phase of extended states in the neigh- 
borhood of the band center. [19] In contrast to the pseudo-random regime, 
it is the average wavelength of the potential oscillations that diverges in the 
thermodynamic limit. Because of the slower divergence of the oscillation wave- 
length, this regime is usually named as weakly aperiodic. The uniform case is 
recovered at z/ = 0. Recently, the effect of this kind of aperiodicity on the 
phonon and magnon ID modes has attracted a renewed interest. [26,27] 

The time-dependent Schrodinger equation governs the dynamics of the Wan- 
nier amplitudes, [41] 



where the Planck constant h is set to unity. 



(3) 



We solve numerically Eq. (3) to study the time evolution of an initially Gaus- 
sian wave-packet centered at site mo. 



'ipmit = 0) = Aexp 



4A2 



(4) 



where A is the normalization constant and we set A = 1 hereafter. Once 
Eq. (3) is solved for the initial condition (4), we compute the projection of 
the mean position of the wave-packet (centroid) and the average velocity v{t) 
along the field direction 



(5) 



as well as the spread of the wave-function (square root of the mean-square 
displacement) 



V m 

where m{t) = Em"^ |V'm.(^)P- 



(6) 



In addition, we use the exact diagonalization of the Hamiltonian (1) in the 
absence of the electric field {U = 0) to obtain the eigenvalues and eigenvectors, 
and to calculate the do conductance G{E) as a function of energy E. For 
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noninteracting electrons and within the hnear response theory at T = OK, 
G{E) is given by the Kubo-Greenwood formula [42] (in units of e^/h) 



2tt 

G{E) = -Y^\{a\xH-Hx\b)\' 

a,b 

X 5{E - Ea)5{E - Eb) . (7) 

Here, the polarization is taken in the direction of the x-axis, Q is the area of 
the system, x is the position operator, and the indices a and b label the eigen- 
states. We also compute the participation number of a normalized eigenstate 
as follows 



P{E) 



(8) 



For extended states, P{E) is expected to scale linearly with the number of 
sites. More specifically, we will focus on the dc conductance and the partici- 
pation number near the band center, averaged over a narrow energy stripe at 
the band center [-Wc/2, Wc/2]: 

E G{E), (9a) 

\E\<Wc/2 

e = ^ E PiE), (9b) 

''^^ \E\<Wa/2 

where A^^; is the number of eigenstates within the stripe. 



3 Results and Discussion 



In all our numerical simulations, we set the parameters V and na of the ape- 
riodic potential toV — 3 and ttq; = 1, while considering the exponents and 
Uy as variables. When calculating G{E), Eq. (7), the (5-functions were replaced 
by 6{E) = {l/AE)Q{AE/2 - \E\), with AE = 0.1, 9 being the Hcaviside 
step-function. The magnitude of the energy stripe in Eqs. (9a) and (9b) was 
chosen Wc — 0.2. 



3.1 Zero-field electronic states 



Firstly, wc discuss the character of the electronic states of the model under 
study. Figure 1(a) displays the averaged conductance G at the band center 
[E = 0) as a function of u = = I'y, obtained after numerical diagonalization 
of the Hamiltonian (1) in the absence of the external electric field {U — 0). To 
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avoid strong fluctuations of tlie conductance related to particular realizations 
of the aperiodic potential, we further averaged the results over a narrow range 
of potential amplitudes {AV = 0.05). As we observe from Fig. 1(a), for a 
slowly varying aperiodic site potential {u^, Uy < 1) the conductance grows on 
increasing the system size. This is a signature of macroscopic conductance 
in the thermodynamic limit. In contrast, the conductance is almost null for a 
pseudo-random potential [vx, > '^)-, in agreement with previous calculations 
for localized eigenmodes. [42] A similar trend is found in the behavior of the 
averaged participation number ^, as shown in Fig. 1(b). For weakly aperiodic 
potentials, ^ grows with the system size, while it becomes size independent 
for strongly aperiodic (pseudo-random) systems. 
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Fig. 1. (a) Averaged dc conductance G, Eq. (9a), and (b) averaged participation 
number ^, Eq.(9b), versus the exponent v = = Vy calculated for different sys- 
tem sizes N . The average was performed considering a narrow range of potential 
amplitudes (Al/ = 0.05) to reduce statistical fluctuations. 
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Fig. 2. Finite size scaling of the normalized participation number ^/N"^ for the states 
in the vicinity of the band center calculated for various values of the exponent v 
(shown in the legend). The inset presents the finite-size scaling of the normalized 
conductance G/N. The different character of the finite-size scaling of both quantities 
on increasing the exponent u signals the change in the nature of electronic states in 
the band center from delocalized {v < 1) to localized ones {v > 1). 

The correlations found in the behavior of the averaged conductance G and the 
participation number ,^ as a function of the exponent u indicate that the model 
under study undergoes a metal-insulator transition associated with a change 
in the nature of the one-electron eigenstates. To give further support to this 
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finding, wc performed a finite-size scaling analysis of these two quantities. In 
Fig. 2 we plotted the size dependence of the normalized participation number 
^/N"^ for various values of the exponent v. As it can be noticed, ^/N"^ does not 
depend of N in the regime of weak aperiodicity < 1), i.e., the participation 
number itself ^ is proportional to the total number of sites N"^. This is a clear 
indication that the eigenstates at the band center are delocalized over a finite 
fraction of the lattice. On the contrary, for strong aperiodic potentials (z/ > 1), 
^/N"^ decreases with the system size, signaling the wave-function localization. 
The size dependence of the normalized conductance G/N m the center of the 
band, shown in the inset of Fig. 2, corroborates the conclusion above. In the 
slowly varying aperiodic regime G (x. N leading to a macroscopic transport 
supported by extended states. Such a behavior is expected for a regular square 
lattice, where the conductance at the band center should be proportional to 
the number of conducting channels, i.e., to the hnear size of the lattice N. In 
the limit of strongly aperiodic lattices, the normalized conductance decreases 
with the system size, thus indicating the absence of macroscopic transport 
in the thermodynamic limit. This is in agreement with the standard scaling 
theory of Anderson localization. [9] 



3.2 Dynamics 

We also study the time evolution of a wave-packet initially localized in the 
center of the lattice, mo = {N/2, N/2), to uncover the character of the elec- 
tronic wave-packet dynamics in the present model system. We start our anal- 
ysis analyzing the wave-packet dynamics in the absence of the external field 
{U = 0). We numerically integrate the wave-equation (3) and calculate the 
wave-packet spread a{t) until it reaches a stationary value as a result of mul- 
tiple reflections of the electron on the lattice boundaries. The extended states 
that appear in this model emerge in the regime at which the potential has 
a diverging wavelength in the thermodynamic limit. These are typically non 
scattered modes thus leading to a ballistic wave-packet spread, a(t) oc t. This 
means that, as long as a{t) < N, the data points calculated for different 
system sizes should collapse into a unique straight fine after re-scahng 
a{t) — > cr(t)/N and t — > t/N. In the localized regime, such a collapse is not 
expected. We further calculated the time-dependent participation second mo- 
ment ^{t) = 1/ J2m\'^mit)\^. In general, ^(t) scales linearly with the number 
of sites whenever the system supports extended states. 

In Fig. 3, we plot the results of our calculations. For a weakly aperiodic site 
potential {u^ = Vy = 0.5), we found the above mentioned collapse of all curves 
into a unique straight line for a{t) < N [see Fig. 3(a)]. This allow us to con- 
clude that the wave-packet motion in this case has a ballistic nature, signaling 
the presence of a phase of delocahzed non scattered states. In the case of a 
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pseudo-random site potential (ux — J^y = ^■5), such a collapse is seen only for 
a{t) <^ [see Fig. 3(b)]: the wave-packet motion is ballistic whenever a{t) is 
smaller than the largest localization length. For long times, a collapse is absent 
[see Fig. 3(b)], indicating the wave-packet locahzation. In Fig. 3(c), we show 
results for the time-dependent participation second moment for i^x — J^y — 0.5 
(solid line) and i^x = ^y = 1-5 (dashed line). For v < 1, the participation 
dynamics show a ballistic behavior until the wave-packet reaches the lattice 
boundaries. For i/ > 1 a slower dynamics was obtained. These results are in 
perfect agreement with the wave-packet spread calculations, indicating the 
presence of non scattered modes in the initial wave-packet only in the regime 
of weakly aperiodic potentials {v < 1). These data corroborate our statement 
about the Anderson transition for the underlined 2D aperiodic model, which 
we claimed in the previous section on the basis of the finite size scaling of the 
conductance G and participation number ^ for the states in the vicinity of the 
band center. 




Fig. 3. The normalized wave-packet spread a{t)/N, Eq. (6), as a function of the 
re-scaled time t/N computed for lattices oi N x N = 200 x 200 to 800 x 800 
sites at zero bias (U = 0). (a) A slowly varying (weakly aperiodic) site poten- 
tial (ux = Uy = V = 0.5). The collapse of all curves into a unique curve until 
a{t) < N indicates the presence of non scattered extended states, (b) A pseu- 
do-random potential {ux = Vy = v = 1.5). A similar collapse is absent, signalling 
localization, (c) Time dependent participation ^(t) number versus time t computed 
for A^ X = 1600 x 1600 and Ux = Uy = f = 0.5 (solid line) and Ux = Uy = u = 1.5 
(dashed line). In agreement with the wave-packet spread calculations, the partici- 
pation number spreads ballistically {^{t) oc i^) for < 1 indicating the presence of 
non scattered modes. 

Now, we turn to the wave-packet dynamics of an electron subjected to a uni- 
form electric field {U 7^ 0). It is well known that in disorder- free systems. 
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Fig. 4. Bias driven dynamics of the centroid R{t), Eq. (5), in a pseudo-random 
potential {vx = = 1-5) with N x N = 200 x 200 sites and bias magnitude 
U = 0.5. No signature of coherent Bloch oscillations is seen in this case. The Fourier 
transform R{u!) of the centroid R{t) is broad, supporting the absence of a typical 
oscillation frequency in R{t). 



96 h (a) (n^^.mJ=(N/2.N/2} 
931 
90 1 
87 1 

84 p- ' U=0.5 

93 [ (h) Jii^tnJ^N/2,N/2) 
C~90| 
8?! 

' U=0.75 




0.4 
0.2 




0.2 



5' 





0.2 



0.1 
0, 



Jd) 




:(e) 




-(f) 

— M+H 1 


w 


1 

-.^ 



0.25 0.5 0.75 1 1.25 1.5 
CO 



Fig. 5. Bias driven dynamics of the centroid R{t) in a slowly varying (weakly ape- 
riodic) potential = Uy = 0.5) for two magnitudes of the bias: (a) U = 0.5 and 
(b) U = 0.75. The initial location of the wave-packet is (no,?7io) = {N/2, N/2). 
(d) and (e) Fourier transforms R{uj) of the centroids depicted on panels (a) and 
(b), respectively. Note that R{io) is peaked at a frequency w ~ [/. (c) Same as in 
panel (a), but for the initial condition (no, mo) = {N/2, N/2 + do) with do = 10. 
(f) Fourier transform of the centroid depicted on panel (c). Note that here R{u) 
exhibits a doTiblet structure. 

a uniform electric field causes dynamic localization of the electron and gives 
rise to an oscillatory motion of the wave-packet, the so-called Bloch oscilla- 
tions. [28,29] The size of the segment over which the electron oscillates and 
the period of the oscillations are estimated from semiclassical arguments to 
he Lp = aW/U and tb — Stt/C/, respectively, [43] where W is the width of 
the Bloch band. 

Firstly, we compute the wave-packet centroid R{t) in a pseudo-random po- 
tential {Px = = 1-5) of size N X N — 200 x 200 with the bias magnitude 
U = 0.5. As deduced from Fig. 4, there is no signature of Bloch oscillations in 
this case. Coherent oscillations, which are present immediately after the initial 
wave-packet is released, are quickly destroyed (not shown). The asymptotic 
behavior of the centroid resembles a stochastic motion around some mean po- 
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Fig. 6. Phase diagrams velocity-vs-position computed for: (a) 
U = 0.5 and the initial focation of the wave-packet (no, mo) 



= My = 1.5, 

(iV/2,iV/2); 

(b) i^x = = 0.5, U = 0.5 and the initial location of the wave-packet 
(no, mo) = {N/2, N/2); (c) the same as in (b) except the initial location of the 
wave-packet (no, mo) = {N/2, N/2 + do) with do = 10- For m > 1 the velocity-vs- 
position diagram exhibits incoherent orbits. For u < 1, the Bloch-likc oscillations 
promotes a coherent dynamics through the phase space. In panel (c), the ampli- 
tude of the oscillations varies faster than in panel (b) due to the doublet structure 
generated by the anisotropic gradient of the local potential. 

sition. The Fourier spectrum R{uj) of the centroid, plotted in the right panel of 
Fig. (4), confirms this claim. The spectrum is rather broad, with no character- 
istic frequencies, indicating that R{t) is similar to a white noise signal. Thus, 
for Vx-i^y > 1 the system shows a behavior similar to the standard Anderson 
model, with no signature of coherent Bloch oscillations. 



In the case of a weakly aperiodic potential {ux = = 0.5), the centroid reveals 
an oscillatory amplitude-modulated pattern, as seen in Fig. 5 (panels a, b, and 
c). Its Fourier transform R{uj) displays a well-defined narrow peak, the location 
of which is slightly deviated from the expected value u = U (panels d and 
e); this small shift is caused by the local contribution to the external bias 
produced by the modulated potential. Indeed, the bias at position (no,mo) is 
given by Ues — {Ssn,m/ + U, Sen,m/ SiTT' + U) . The local contribution shall be 
relevant whenever the potential gradient (in appropriate units) is of the order 
of the external bias U . Therefore, the frequency obtained in not exactly uo = U 
[Figs. 5(a) and 5(b)]. The two peak structure depicted in Fig. 5(c) has its origin 
in the distinct contributions given by the anisotropic gradient of the local 
potential to the effective local bias. Once the initial wave-packet is located at 
a non- symmetrical position of the lattice {tiq, itlq) = {N/2, N/2-\-dQ) with do = 
10, the potential derivatives along the orthogonal lattice directions give raise 
to distinct frequency shifts when the electron is forced by a diagonal external 
field. We stress that such frequency sphtting is not found in ID aperiodic 
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Fig. 7. The dynamics of the centroid R{t) and its Fourier transform (left and right 
panel respectively) for = 0.5 and (a) Vy = 1.5 and (b) Vy = 0.7. For both 
and Uy below 1, Bloch-oscillation with a frequency close to the one predicted by the 
semi-classical approach was obtained. However, in the case of mixed pseudo-ran- 
dom and weakly aperiodic components, the pseudo-random character for i/y > 1 
predominates leading to the Anderson localization of the one-electron eigenstates 
and the consequent incoherent dynamics of the wave-packet centroid. 

potentials, being thus a specific feature of higher dimensional systems. 

In Fig. 6, we present the phase space plots, namely velocity versus position . 
All calculations were performed using NxN = 200 x 200. Three representative 
cases are illustrated: (a) u^. = Uy = 1.5, U = 0.5 with the initial location of the 
wave-packet (no, mo) = {N/2, N/2); (b) i^x = = 0.5, U — 0.5 with the initial 
location of the wave-packet (no, mo) = {N/2, N/2); (c) the same as in (b) 
except by initial location of the wave-packet (no, mo) = {N/2, N/2 + do) with 
do = 10. For all cases, the integration time used was t = 1000. We can see in 
(a) that the velocity versus position diagram exhibits orbits with no coherency, 
as it was observed in the centroid dynamics (see Fig. 4). In the other cases, 
the diagrams show a coherent dynamics through the phase space. In the case 
of bimodal Bloch oscillations [see Fig. 6(c)], the wave-packet orbits in phase 
space show a pronounced breathing pattern. This breathing reflects the fact 
that both centroid and velocity display amplitude- modulated envelopes. The 
frequency of this amplitude- modulated envelope depends on the difference 
between the two frequency peaks appearing the Bloch oscillations spectrum 
[Fig. 5(c)]. 

Before concluding, some words concerning the electronic dynamics in systems 
with distinct degrees of aperiodicity along the orthogonal lattice directions 
J^x 7^ In Fig-7, we show the centroid R{t) and its Fourier transform (left 
and right panel respectively) for an initial Gaussian wave-packet located at 
the center of the lattice with aperiodicity exponents = 0.5 and (a) Uy = 1.5 
and (b) Uy = 0.7. For both and Uy below 1, we can see oscillations with 
a frequency close to the Bloch-oscillation frequency predicted by the semi- 
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classical approach. However, in the case of mixed pscTido-random and weakly 
aperiodic components, the pseudo-random character for Uy > 1 predominates 
leading to the Anderson localization of the one-electron eigenstates and the 
consequent incoherent dynamics of the wave-packet centroid. 



4 SummEiry and concluding remarks 

In this work, we studied numerically the dynamics of a single electron wave- 
packet moving in a 2D square lattice with an aperiodic site potential = 
V cos (Tram^^) cos (Tram^^) in the presence of an external uniform electric field. 
For fixed parameters V and na, the exponents and Uy allow to control the 
degree of aperiodicity. 

We have discussed the character of the electronic states of the model in the 
absence of an external electric field. Diagonalizing numerically the Hamilto- 
nian at zero electric field, we computed the dc conductance and the partici- 
pation number for various systems sizes. After a finite-size scaling analysis of 
data close to the band center, we found that, for a weakly aperiodic potential 
{ur^, Uy < 1), both quantities grow on increasing the system size. The participa- 
tion number is proportional to iV^ while the conductance scales proportional 
to A^. Such size dependencies are consistent with the presence of extended 
states and macroscopic transport in the thermodynamic limit. In contrast, 
both participation number and conductance reach a plateau upon increasing 
the size for a pseudo-random potential {I'x-i^y > 1), which is the typical be- 
havior predicted by the usual scaling theory of Anderson localization. Solving 
numerically the time-dependent Schrodinger equation, we found that an ini- 
tially localized wave-packet reveals a ballistic spreading associated with the 
existence of a phase of extended non scattered states in the weakly aperiodic 
limit {i^x-i^y < 1). In contrast, the spread of the wave-packet is bounded for 
pseudo-random realization of the site potential {ux.i'y > 1). These results 
indicate that the underlined 2D aperiodic model under study undergoes a 
metal-insulator transition. 

We have also investigated the interplay between the delocalization effect, pre- 
served by the weakly aperiodic structure, and the dynamic localization, caused 
by an electric field acting on the system. We demonstrated that in weakly 
aperiodic potentials {i^x-i^y < 1); the applied uniform electric field promotes 
sustained Bloch oscillations of the electronic wave-packet. The frequency of 
the oscillations was shown to agree with the prediction of a semi-classical 
approach. [43] In contrast with the Bloch oscillations in ID, the frequency 
spectrum of the centroid time-evolution can exhibit a two-mode structure 
depending on the initial position of the wave-packet. Such two- frequency dy- 
namics reflects the anisotropy of the potential landscape along the principal 
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lattice directions when the electron is forced to oscillate by a diagonal field. 
In phase-space, this phenomenon is signaled by quickly breathing coherent 
orbits. 

5 Acknowledgments 

Work at Alagoas was supported by CNPq-Redc Nanobioestruturas, CAPES 
(Brazihan research agencies) and FAPEAL (Alagoas State agency). Work 
at Madrid was supported by MEC (Project MOSAICO) and BSCH-UCM 
(Project PR34/07-15916). 

References 

[I] G. Khitrova, H. M. Gibbs, F. Jahnke, M. Kira, and S. W. Koch, Rev. Mod. 
Phys. 71, 1591 (1999). 

[2] See contributions to Optics of Quantum Dots and Wires, eds. G. W. Briant and 
G. S. Solomon (Boston: Artech House, 2005). 

[3] T. Brandes, Phys. Rep. 408, 315 (2005). 

[4] J. R. Tischler, M. S. Bradly, Q. Zhang, T. Atay, A. Nurmiko, and V. Bulovic, 
Org. Electronics 8, 94 (2007). 

[5] E. L. Albuquerque and M. G. Cottam, Phys. Rep. 376 225 (2003); Polaritons 
in Periodic and Quasiperiodic Structures (Elsevier: Amsterdam, 2004). 

[6] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[7] N. F. Mott, J. Non-Cryst. Solids 1, 1 (1968). 

[8] N. Mott and W. D. Twose, Adv. Phys. 10, 107 (1961). 

[9] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, 
Phys. Rev. Lett. 42, 673 (1979). 

[10] J. C. Flores, J. Phys.: Condens. Matter 1, 8471 (1989). 

[II] D. H. Dunlap, H.-L. Wu, and P. W. Phillips, Phys. Rev. Lett. 65, 88 (1990). 
[12] P. W. Phillips and H.-L. Wu, Science 252, 1805 (1991). 

[13] A. Sanchez and F. Dommguez-Adame, J. Phys A:Math. Gen. 27, 3725 (1994); 
A. Sanchez, E. Macia, and F. Dommguez-Adame, Phys. Rev. B 49, 147 (1994). 

[14] E. Diez, A. Sanchez, F. Dommguez-Adame, Phys Rev. B 50, 14359 (1994); 
F. Dommguez-Adame, E.Diez, and A. Sanchez, Phys. Rev. B 51, 8115 (1995). 



14 



[15] F. A. B. F. de Moura and M. L. Lyra, Phys. Rev. Lett. 81, 3735 (1998); 84, 
199 (2000); Physica A 266, 465 (1999). 

[16] F. M. Izrailev and A. A. Krokhin, Phys. Rev. Lett. 82, 4062 (1999). 

[17] V. Bellani, E. Diez, R. Hey, L. Toni, L. Tarricone, G. B. Parravicini, F. 
Dommguez-Adame, and R. Gomez-Alcala, Phys. Rev. Lett. 82, 2159 (1999). 

[18] U. Kuhl, F. M. Izrailev, A. A. Krokhin, and H.-J. Stockmann, Appl. Phys. Lett. 
77, 633 (2000). 

[19] S. Das Sarma, S. He, and X. C. Xie, Phys. Rev. Lett. 61, 2144 (1988); Phys. 
Rev. B 41, 5544 (1990). 

[20] D. R. Grempel, S. Fishman, and R. E. Prange, Phys. Rev. Lett. 49, 833 (1982). 

[21] M. Griniasty and S. Fishman, Phys. Rev. Lett. 60, 1334 (1988). 

[22] D. J. Thouless, Phys. Rev. Lett. 61, 2141 (1988). 

[23] H. Yamada, Phys. Rev. B 69, 014205 (2004). 

[24] E. Macia, Rep. Prog. Phys. 69, 397 (2006). 

[25] E. Macia and F. Dommguez-Adame, Electrons, Phonons and Excitons in Low 
Dimensional Aperiodic Systems, Editorial Complutense, Madrid (2000). 

[26] F. A. B. F. de Moura, L. P. Viana, and A. C. Frery, Phys. Rev. B 73, 212302 
(2006). 

[27] F. A. B. F. de Moura, Eur. Phys. J. B 58, 389 (2007). 
[28] F. Bloch, Z. Phys. 52, 555 (1927). 

[29] C. Zener, Proc. R. Soc. London, Scr. A 145, 523 (1934). 
[30] L. Esaki and R. Tsu, IBM J. Res. Div. 14, 61 (1970). 

[31] G. H. Wannier, Phys. Rev. 100, 1227 (1955); Phys. Rev. 101, 1835 (1956); 
Phys. Rev. 117, 432 (1960); Rev. Mod. Phys. 34, 645 (1962). 

[32] D. H. Dunlap and V. M. Kenkre, Phys. Rev. B 34, 3625 (1986). 

[33] K. Leo, P. Hairing, F. Briiggemann, R. Schwedler and K. Kohler, Solid State 
Commun. 84, 943 (1992). 

[34] K. Leo, Semicond. Sci. Technol. 13, 249 (1998). 

[35] V. Agarwal, J. A. delRio, G. Malpucch, M. Zamfirescu, A. Kavokin, D. 
Coquillat, D. Scalbert, M. Vladimirova and B. Gil, Phys. Rev. Lett. 92, 097401 
(2004). 

[36] F. Domfnguez-Adame, V. A. Malyshev, F. A. B. F. de Moura, and M. L. Lyra, 
Phys. Rev. Lett. 91, 197402 (2003). 



15 



[37] F. A. B. F. dc Moura, M. L. Lyra, F. Dommguez-Adame, and V. A. Malyshev, 
J. Phys.: Condens. Matter 19, 056204 (2007). 

[38] F. Dommguez-Adame, V. A. Malyshev, F. A. B. F. de Moura, and M. L. Lyra, 
Phys. Rev. B 71, 104303 (2005). 

[39] E. Diaz, F. Dommguez-Adame, Yu. A. Kosevich, and V. A. Malyshev, Phys. 
Rev. B 73, 174210 (2006). 

[40] H. N. Nazareno and P. E. de Brito, Phys. Rev. B 60, 4629 (1999). 

[41] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56 1469, (1993). 

[42] M. Tit and M. Schrcibcr, J. Phys.: Condens. Matttcr. 7, 5549 (1995). 

[43] N. W. Ashcroft and N. D. Mcrmin, Solid State Physics (Saunders Colege 
Publishers, New York, 1976), P. 213. 



16 



